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A theory for boson-assisted tunneling via randomly distributed resonant states in a layered metals 
is developed. As particular examples, we consider the electron-phonon interaction and the interac- 
tion between localized and conduction electrons. The theory is applied to explain a non-monotonic 
variation of the out-plane resistivity with temperature observed in quasi-two-dimensional metals. 
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I. INTRODUCTION 



Electron transport in layered metals exhibits a 
qualitatively different behavior of the in-plane (p a b) 
and out-of-plane (p c ) resistivities: whereas the tem- 
perature dependence of p a b is metallic-like, that of 
p c is either insulating-like or even non-monotonic. 
This behavior is observed in various materials, 
such as high-temperature superconductors^, sodium 
and bismuth cobaltate oxides^^, the layered per- 
ovskite S^RuO^^, dichalogenidesS, graphite^, or- 
ganic metals^, and other systems. It is quite remarkable 
that the c-axis resistivity behaves similarly in materials 
with otherwise drastically different properties, ranging 
from weakly (graphite) to strongly (S^RuCU) correlated 
Fermi liquids and then to non- Fermi liquids (HTC, cobal- 
tate oxides), and with Fermi energies spanning the inter- 
val from a few eV (in most layered metals) to 25 meV 
(in graphite). Also, the magnitude of p c varies from a 
few mf2-cm (e.g., in Si^RuO^) to a few fi-cm (e.g., in 
organics and (Bii^Pbr^S^C^Os). The great variety 
of systems and associated energy scales suggests that the 
origin of this effect is not related to the specific properties 
of any of the compounds but rather to what they have 
in common, i.e., strong anisotropy. 

At the level of non- interacting electrons, layered sys- 
tems are metals with strongly anisotropic Fermi surfaces. 
A simple but instructive model is that of free motion 
along the layers and nearest-neighbor hopping between 
the layers. In this model, the single-particle spectrum is 



e k = kj|/2m afc + 2J C (1 - cosk z d) 



(1) 



where k|| and k z are in the in-plane and c-axis compo- 
nents of the momentum, respectively, m a b is the in-plane 
mass, J c is the hopping matrix element in the c-axis direc- 
tion, and d is the lattice constant in the same direction. 
Strong mass anisotropy is guaranteed by the condition 
m a b <C m c = 1/2 J c d 2 .(We set fi = fcs = 1 through the 
rest of the paper.) For Ep < 2J C the Fermi surface is 
closed (as it is in graphite); for Ep > 2J C the Fermi sur- 
face is open (as it is in the majority of layered materials). 
For weakly coupled layers (Ep ^s> J c ), the equipotential 
surfaces are "corrugated cylinders" with slight modula- 
tion along the c-axis (see Fig.l). 

Transport in metals is commonly described via the 




FIG. 1: Fermi surface corresponding to Eq.([TJ with Fermi 
velocity vectors at two different points. 



Boltzmann equation for the distribution function of elec- 
trons (/) in the phase space. To study the linear response 
conductivity, it is sufficient to consider a weak electric 
field. In this case, / = /o + /i, where the non-equilibrium 
correction f\ to the Fermi function fo satisfies 



d h i n f 
^+v.V r/l 



7[/ 1 ] = -^evE. (2) 



Here J is a linearized collision integral and v = Vk£k- 
The tensor of the dc electrical conductivity is given by 



d 3 p ( df Q 



(2<T 



^jv^pK-^p), (3) 



where the integration goes over the Brillouin zone. For 
all known types of the inelastic interaction (electron- 
phonon, electron-electron, electron-magnon, etc.), the 
operator I^ 1 decreases with the temperature; thus all 
components of the conductivity tensor must decrease 
with T as well (although the particular forms of the T 
dependence may be different for different components). 
This is not what the experiment shows. 

The situation when Bolztmann equation fails so dra- 
matically is rather unusual, given that other examples 
of its breakdown, e.g., weak localization and Altshuler- 
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Aronov effects in the conductivity, indicate some non- 
trivial quantum interference effects. The fact that this 
breakdown occurs only in the c-direction, makes the sit- 
uation even more puzzling. A vast amount of litera- 
ture, addressing various aspects of this problem, is ac- 
cumulated at this point. The proposed models can be 
roughly divided in two groups. The first group is trying 
to find the explanation of the anomalous c-axis trans- 
port within a model of an anisotropic metal in which 
an electron interacts with potential impurities, phonons, 
and/or other electrons. The breakdown of the Boltzmann 
equation is associated either with Anderson localization 
in the c (but not in the in-plane) direction or from the 
"coherent-incoherent crossover" , which is believed to oc- 
cur when the tunneling time between the layers {J^ 1 ) 
becomes longer than the inelastic time (ri n )i2 or the ther- 
mal time (T -1 )^ 3 -, i.e., when J c min{Tj n , T -1 } < 1. How- 
ever, as far as Anderson localization is concerned, it was 
shown that it occurs only for exponentially small J c [J c < 
r^ 1 exp {—Eft c \) , where r e i the elastic scattering time] 
and then only simultaneously in all direction o 14 ^ 15 ' 16 . 
The breakdown of the Boltzmann equation due to the 
coherent-incoherent crossover for the electron-phonon in- 
teraction has been reconsidered in our recent paper—, 
where we have shown that the Prange-Kadanoff deriva- 
tion of the Boltzmann equation 1 ^ applies equally well 
to the anisotropic case. The only condition for the va- 
lidity of the Boltzmann equation for both isotropic and 
anisotropic cases is the large value of the "dimensionless 
conductance" EpTi n , where T- m is the inelastic scatter- 
ing time, regardless of whether J c min{Ti n , T -1 } is large 
or small, provided that the Migdal parameter s/vp (s is 
the sound velocity) is small over the entire Fermi surface, 
which is the case for Fermi surfaces of the type shown in 
Fig. [TJ This argument can be readily extended to any sit- 
uation when the self-energy is local, i.e., independent of 
the electron's momentum, and applies not only to Fermi 
liquids but also to non- Fermi liquids. Therefore, for a 
large class of non-perturbative interactions the coherent- 
incoherent crossover does not occur, and the c-axis re- 
sistivity is supposed to maintain its metallic character. 
On the experimental side, the existence of the coherent- 
incoherent crossover has been questioned by a recent ob- 
servation of angular magneto resistance oscillations in a 
layered organic metat^ well above the temperature of the 
expected crossover. Moreover, the coherent-incoherent 
crossover on its own does not account for a non-metallic 
conductivity, as in this scenario the conductivities in both 
the coherent and incoherent regimes are proportional to 
the same scattering time2£j2I and thus exhibit a metallic 
T dependence. 

Along a similar line of reasoning, it was suggested 
that the non-metallic c-axis transport is related to the 
Fermi-to-non-Fermi-liquid crossover (understood as a 
smearing of the quasiparticle peak in the spectral func- 
tion) , which occurs at high enough temperatures in some 
material o 5 ' 22 . Indeed, the angle-resolved photoemission 
experiments on sodium and bismuth cobaltate oxides 



show that the temperature at which the quasiparticle 
peak is smeared is of the same order as temperature Tm 
at which p c exhibits a maximum 5 -^. Although this coin- 
cidence is suggestive, it runs against the Prange-Kadanoff 
argument 1 ^ which shows that the existence of quasipar- 
ticles is not a pre-requisite for Boltzmann- like transport. 
Also, it is not clear at the moment whether the rela- 
tion between smearing of the quasi-particle peak and the 
maximum in p c is common for all materials. For example, 
Tm in the HOPG graphite is low enough (~ 40 K)i£, so 
that quasiparticles are still well-defined at T ~ Tm- Fi- 
nally, the model of Fermi-to-non-Fermi-liquid crossover 
on its own does not account for the anomalous trans- 
port, especially given that the in-plane resistivity shows 
no dramatic signs of this crossover. 

Within the first group is also the polaron mode l 2526 , 
which is capable of a quantitative description of the ex- 
periment, given that polarons are stable. The latter as- 
sumption, however, requires a very strong anisotropy of 
the phonon spectrum, which needs to be more anisotropic 
than the electron one, which is not the case at least in 
some representative materials, e.g., in Si'2Ru04 2 — . Fi- 
nally, there is an explanation based on the zero-bias 
anomaly in tunneling between the layers, resulting from 
a suppression of the single-particle density of states of 
in-plane electrons, e.g., via a pseudogap mechanism in 
high T c cuprates2Sj2£. However, the zero-bias mecha- 
nism is only valid on the incoherent side of the coherent- 
incoherent crossover which, as we argued earlier, does not 
occur for a large class of interaction. 

The second group assumes that the primary reason of 
the anomalous behavior of the c-axis resistivity is trans- 
port through inter-plane defects. The temperature de- 
pendence is introduced either phenomenologicall y 29 ) 30 or 
through thermal occupation of electron states in the con- 
ducting layer o 9 i 31 ' 32 , or else, as it was done in the short 
version of this paper—, through a phonon-assisted tun- 
neling. Whether inter-plane disorder is the reason for 
the anomalous c-axis transport in all layered materials is 
not clear at the moment. Currently, there is a number of 
observations pointing at the role of inter-plane disorder 
as a mediator of c-axis transport. For example, recent 
experiment 3 - has shown that the radiation damage of an 
organic metal reduces rather then increases p c at tem- 
peratures around Tm- Also, variability of p c in graphite 
samples prepared in different ways (including the natu- 
ral ones) provides an indirect argument for the role of 
inter-plane disorder. It seems worthwhile to explore the 
consequences of a model involving inter-plane disorder 
and this what we will do in this paper. 

We assume that, in addition to normal impurities, a 
layered crystal also contains a number of resonant impuri- 
ties located in the inter-planar space. An electron moves 
between the layers via two mechanisms: the first one is di- 
rect tunneling, augmented by scattering at normal (non- 
resonant) impurities, phonons, etc., and the second one is 
resonant tunneling through defect sites. Direct tunneling 
between the layers forms a band state smeared by vari- 
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ous scattering processes, which include the non-resonant 
part of scattering by inter-plane defects. Transport of 
these states is described by the Boltzmann conductivity 
erf which has a metallic temperature dependence. As 
resonant tunneling opens a new channel of conduction, 
the total conductivity can be described by a phenomeno- 
logical formula 7 ' 30 ! 34 

a c = a? + af«, (4) 

where a T c es is the resonant-impurity contribution. On the 
other hand, the in-plane conductivity remains largely un- 
affected by inter-plane disorder, so that a ao — cr^ b . Re- 
gardless of a particular tunneling mechanism, the reso- 
nant part al cs increases with T. Consequently, the band 
channel, which is weak to begin with due to a small value 
of the inter-plane transfer J c , is short-circuited by the res- 
onant one at high enough temperatures. Accordingly, a c 
goes through a minimum at a certain temperature (and 
p c = c^ 1 goes through a maximum). More generally, p c 
may exhibit a variety of behaviors , discussed in Sec lIVI 
In this paper, we develop a microscopic theory of res- 
onant tunneling through a wide band of energy levels 
positioned randomly in between two conducting layers 
and coupled to a fluctuating field of bosonic excitations 
fSec lII Aj ). Two particular examples of such a field, con- 
sidered here, are phonons and the dynamic Coulomb field 
of all conduction electrons in a crystal. The case of 
phonon-assisted tunneling through a single junction was 
considered before in Refs. [35ll36| | Our formalism repro- 
duces the general results of Refs. (3f|;36], a s well as the 
low— T behavior of the conductivity, found in Ref. [35| . 
In Sec. IIIB1 we obtain a detailed, non-perturbative ex- 
pression for a r c es for the electron-phonon interaction and 
showed that a r c cs saturates at temperatures higher than 

r s = Xuj d , (5) 

where wb is the Debye frequency and A is the dimcn- 
sionless coupling constant for on-site electrons. This 
prediction is important for discriminating the phonon- 
assisted mechanism against other effects. A strong on- 
site Coulomb interaction was shown to result in a Kondo 
anomaly in the tunneling conductanc e 37 ' 38 . The effect 
of the Coulomb interaction between on-site and conduc- 
tion electrons was considered by Matveev and Larkin 39 
in the context of a threshold singularity in the non-linear 
current- voltage characteristic, observed in Ref. [4p| . The 
new element of this work is that we consider the effect 
of the Coulomb interaction between on-site and conduc- 
tion electrons on the linear tunneling conductance both 
in the ballistic and diffusive regimes of conduction elec- 
trons' motion fsec. Ill Cj) . The most interesting result of 
this analysis is the scaling behavior a r c es oc T VB in the 
ballistic regime, where t\b is the dimensionless coupling 
constant for the Coulomb interaction. In Sec lIIIl we ana- 
lyze the results. In Section [TVl we compare our theory to 
the experiment and show that equation ^ describes well 



the non-monotonic T-dependence of p c in S^RuO^and 
k-(BEDT-TTF) 2 Cu(SCN)25 4 . Our conclusions are given 
in SecEl 

II. BOSON- ASSISTED TUNNELING 
A. General formalism 

We consider tunneling through a single resonant impu- 
rity located in between two metallic layers. To account 
for a finite concentration of such impurities, we will av- 
erage the result with respect to the positions of the res- 
onant center in real space and also within the energy 
band. In doing so, we neglect the effects of interference 
between different resonant centers, as well as between res- 
onant and non-resonant scattering. Tunneling through 
more than one impurity was considered theoretically in 
[4l]and observed experimentally in thick (> 1 nm) tunnel- 
ing junctions 4 ^, but it is less likely to occur in tunneling 
through thin ( < 1 nm) inter-plane spacings, so we will 
disregard such a possibility here. 

We use the tunneling Hamiltonian description 

H = H a + H d + H c + H ac + H C( i + -ffinoi • (6) 
The free part of the Hamiltonian 

k ll k l 

represents metallic layers. The tunneling part of the 
Hamiltonian 

H ac =535ca(k||)(4| | c + c t a k | l ), (8) 

k 

ffcd = $>cd(k| | )(4„c + (9) 
k n 

describes hopping on and off the resonant site, located at 
point i"i = (r|| = 0,Zi) in between the layers. Dynamics 
of the resonant level is accounted by 

H c = (eo + (j>{?i,t))c ] c, 

where the time-dependent operator <j){r, t) describes the 
fluctuations of the electrostatic potential. Although a 
general result can be obtained for an arbitrary (but Gaus- 
sian) field (f>, we will be mostly interested in the case when 
this field arises due to fluctuations in the positions of ions 
and electrons 

4>(r,t) =0P h (r,i)+0 c (r,i). (10) 

The contribution arising from the crystal degrees of free- 
dom is the potential produced by the displacement wave 

0P h (r, t) = J2 a q (^e l ^ t - iq - r - b ^ t+ ^- r ) , (n) 
q 



4 



where 6 q is the phonon creation operator, a q is the vertex 
of the electron-phonon interaction, and u> q is the phonon 
dispersion. The electronic part of the potential is ex- 
pressed through fluctuation of the electron density p in 
the entire crystal, including the two layers accounted for 
in H a and Hd 

$ e (r,t) = J dr'V (\r-r'\)p(r',t), (12) 

where V (\r\) = Aire 2 /\v\. 

We assume that the potential localizing the electron is 
of a short range and therefore, when the site is empty, it 
does not affect the motion of conduction electrons. This 
situation is different from that considered in the context 
of the Fermi edge singularity in resonant tunnelin g 39 ' 46 , 
where the Coulomb field of the empty, charged site scat- 
ters conduction electrons. We also assume that the on- 
site Coulomb repulsion is so large that the double occu- 
pancy of the resonant center is forbidden, and the spin 
degree of freedom does not play any role in tunneling. 
The Kondo effect in tunneling through a resonant impu- 
rity was considered in Refs. [37[. 

The last term in the Hamiltonian, H lnc \, consists of the 
frec-phonon part 



neling. By analogy the effect of the electron-electron 
interaction can be called electron-assisted tunneling. We 
show that the problem of assisted tunneling allows for 
an exact solution for any bosonic field <f>, given its fluc- 
tuations in time are Gaussian. Remarkably, the effect of 
different bosonic degrees of freedom (such as phonons, 
plasmons, or diffusive density fluctuations), affecting the 
tunneling probability at various temperatures, can be in- 
corporated into a single formula. 

The tunneling current is obtained from the balance 
equation 



l=eT W k||)k{ <(l 



k||,k' 



Wk {| ,k„<(l 



•k' ) , 

where Wkn.p,, is a transition rate through the resonant 



center and 



L/R 



are the distribution functions for elec- 



trons in the left and right layers, respectively. In the 
linear response regime, the conductance of a bi-layer tun- 
neling junction is 



G 



dede'W e 



dn e < 
n € >) + -z-rn t 
oe' 



(15) 



H 



P h 



w q & q & q • 



(13) 



and of the part describing the dynamics of all mobile 
electrons in the crystal. This dynamics controls the be- 
havior of the density fluctuations, p(r, t) . 

The effect of the electron-phonon interaction on tun- 
neling is traditionally referred to as phonon-assisted tun- 



where e = eit| j , e' = e^' and is the dispersion for the 
in-plane motion. The transition rate can be expressed in 
terms of the resonant tunneling amplitude 3 ^ 



e , = Jim hul e ,{t)U e>e ,(t)), 



where the transition amplitude is given by 



(16) 



U e ,e>(t) = -<7 C aff cd e i(£ - £,)t / dne'^ ^'-^ r dT 2 e i ^- e - ir ^Texp\-i P $(t)dt\ 

Jo JO I Jt 2 ) 

Integrating out the Gaussian field </>, one gets 

/•OO /"OO 

W e / = g 2 ca g 2 cd / dt ie ~^-^ / Ma^- 4 ^-^ 3 ^ -^-^!,^) , 



(17) 



(18) 



where T = T L +T R , T L = V g% a 5{e - e P|1 ), T R = V g 2 d S(e - e Pn ), 



V(ti,t 2 ,t 3 ) = exp ji J (u),Zi) 



iuita — t 2 ) — e iujta — e~ iuJt2 — e iutl + e iuJ ^ tl+t3 ' 1 + e M*i-*2) _ e <w(t3-t 2 +ti) 



coth ^— j ^e iuJt3 + e~ iuJt2 — 2 + e iLOtl — e <w (*i+*3) _ e i«(ti-t2) _|_ e i^(t 1 +t 3 -t 2 )^j 



(19) 



Here S (uj, Zj) is the local spectral function of potential fluctuations, as measured at the resonant site 



S (w,Zi) 



it J (2tt) 2 



ImD fl (w,q||, (20) 
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where the retarded propagator of <j> is 



D R (w,q lb z,z')=-i J dt I d 2 r ll e i ^ t -^n)( j, ( r|| , z , t ) (0, z', 0) ) 



(21) 



r 



As we assume that the metal is translationally invariant 
in the in-plane direction but periodic in the z-direction, 
the spectral function, in general, depends on Zj. In the 
absence of interactions, Eq. (fT5|) reproduces the Breight- 
Wigner formula. For the electron-phonon interaction, <j) 
is a deformation potential at the resonant site. In that 
case Eq. (fl"8|) reproduces the result of Ref. 

From now on, we consider the case of a resonant- 
impurity band, assuming that the resonant centers are 
randomly distributed over the inter-layer spacing while 



their energies are uniformly distributed in the interval 
Ef, around the Fermi energy. Assuming that 3> T 
and averaging over and £o, we simplify Eq.(TT5]) further 
to 



G- 



where 



W(e) 



1-coth ( — )+ — k-. r 

V2r^2Tsinh 2 (^) 



(de), (22) 



J 



W(e) = Wo / dtexp 



iet - 



duj 



(1 — cos(wi)) coth ( — i sin(wi) 



(23) 



r 



Here Wq is the resonant transition probability in the ab- 
sence of interaction and S m (u>) = S (u>,d/2) is the spec- 
tral function at the position of the most efficient resonant 
center, i.e., in the middle of the spacing between the lay- 
ers. The local spectral function S m {oj) contains all the 
information about the fluctuating field. 

Eqs. (|22l23[) describe the average conductance for a sin- 



gle resonant impurity. Since a layered metal can be 
viewed as sequence of tunneling junctions connected in 
series, the resonant tunneling conductivity of the whole 
crystal is related to the bi-layer conductance by the 
Ohm's law: er res = Gd. Performing integration over en- 
ergy in Eq. (|2"2"|) . we arrive at the general result for the 
boson-assisted tunneling conductivity: 



J 



= 0"cl 



DC 

dt- 



sinlr (nTt + iO) 



exp 



du 

— b m (uj) 



(1 — cos(cji)) coth 



(-) 

\2TJ 



— i sin(wi) 



(24) 



r 



where a c \ ~ 7re Tni mp aod/ Eir* is a resonant conductiv- 



ity of free electrons, 



''imp 



is the number of resonant im- 



purities per unit volume, and ao is the localization radius 
of a resonant state. 



We note a similarity between our result (|24p and those 
of two other problems: the zero bias anomaly in the tun- 
neling current^ and "dissipative localization" of a par- 
ticle coupled to the thermal bath [the Caldeira-Leggett 
(CL) model^]. Indeed, all of these problems are re- 
lated and describe different variations of the Mossbauer 
effect 44 . In the case of tunneling through a resonant im- 
purity, its occupied and empty states play the role of two 
pseudo-spin states of the CL particle, whereas the fluctu- 



ating electrostatic potential is analogous to the thermal 
bath of harmonic oscillators. While this similarity is not 
formally obvious for tunneling through a single resonant 
center, it becomes clear after averaging over the ensemble 
of levels. Despite the similarity, there are also some dif- 
ference with the CL model. In the CL model, the effect of 
the environment is incorporated phenomenologically, via 
the parameters of noise correlation function. A more mi- 
croscopic resonant-tunneling model treats explicitly fluc- 
tuations of the electrostatic potential, incorporated in the 
spectral function S m (co). As the temperature changes, so 
does the characteristic frequency scale and, consequently, 
the most efficient bosonic mode. This results in a num- 
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ber of intermediate regimes of assisted tunneling, each of 
them corresponding to its own effective CL model. 

The next two Sections are devoted to a quantitative 
analysis of some of the most common sources of interac- 
tion: the electron-phonon and electron-electron ones. 

B. Phonon-assisted tunneling 

We assume that phonons are isotropic and that 
the electron-phonon interaction is of the deformation- 
potential type. Although phonon modes in real lay- 
ered metals are anisotropic, this anisotropy is still 
much weaker than the anisotropy of electron spectra. 
Anisotropy of phonon modes can be incorporated into 
the theory without any difficulties. The retarded corre- 
lation function of the deformation potential is given by 

D R {u, q )=al— ^ - , (25) 

I 



where a 2 = A 2 g 2 /pw q , A is the deformation potential 
constant, and p is the atomic mass density. In the elastic 
continuum model, the spectral function is translation- 
ally invariant in the z-direction. For acoustic phonons 
w q = sq6 (q D - q) , Eg. (1201251) ) gives 

S^ h (cu) = -\d(u JD -^4- > (26) 

uj d 

where 

A = A 2 c4/psV (27) 

is the dimensionless coupling constant for localized elec- 
trons, u>d = sqjj, and 9 (x) is the step-function. The 
interaction with acoustic phonons corresponds to the 
super-Ohmic regime of the CL model. Using Eqs. 
f)24|26[) . we obtain 



-<j c \ I dt 



iirTH 



/(*) = 



WD 

dw — 



D 



sinh 2 (nTt + iO) 
((l-cos(wt)) 



~ A/(t) , where 



coth ( — — 
2T. 



(28) 
(29) 



In the absence of the electron-phonon interaction (A = 0), 

(Trcs 0"el • 

Notice that the electron-phonon interaction is much 
stronger for localized electrons than for thewe band ones. 
Indeed, the dimensionless coupling constant for bulk elec- 
trons which determines, e.g., the mass-renormalization, is 
of order of unity for most metals: ( = a 2 qD (gz>) qo/vFS ~ 
1. The coupling constant between localized electrons and 
phonons is larger than £ by at least the Migdal param- 
eter: A ~ ( (kpd) (vp/s) ^> 1. Therefore, one needs to 
consider a non-perturbative regime of phonon-assisted 
tunneling. Analyzing Eg. (f2"B")) (for details, see Appendix 
|A"]) , one finds that resonant tunneling is exponentially 
suppressed at T = 0: 



cr rcs (T = 0) = cr ol e 



-A/2 



(30) 



At finite T, we find 

f£e. = | exp(-f + |(£) 2 ) ,^«T«^ 



(Tel 



{ l-ffi , T > Xu) D . 



(31) 

The low-temperature limit (T <C W/j/VA), reproduces 
the known result of Ref. [35j. As we see, er res increases 



with T, resembling the zero-bias anomaly in disordered 
metals and Mossbauer effect. At high temperatures 
(T ^> T s = \lod), cr r es saturates at the non- interacting 
value ((Tei), because thermal activation cannot make the 
conductivity larger than in the absence of phonons. No- 
tice that, in contrast to the phenomenological model of 
Ref. [3(3], there is no simple relation between the In- 
dependences of erf and cr rcs . 

The temperature dependence of the tunneling conduc- 
tivity also resembles the polaronic behavio r 25 i 26 i 45 . In- 
deed, phonon-assisted tunneling can be interpreted in 
terms of a local polaron formation. However, we empha- 
size that polarons are stable here because the velocity 
of on-site electrons is determined by the tunneling width 
and thus very small. In the absence of resonant impu- 
rities, a metal with a large Migdal parameter vf/s ^> 1 
cannot sustain stable polarons, as they emit phonons and 
decay. 



C. Electron-assisted tunneling 

The effect of the interaction between on-site and con- 
duction electrons on resonant tunneling is qualitatively 
similar to the effect of phonons, as a tunneling electron 
drags a surrounding cloud of electrons. This mechanism 
is also known as the "electronic polaron"—. 

The on-site Coulomb potential is produced by all mo- 
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bile electrons in a metal. The Green's function of elec- 
trons on a lattice is not translationally invariant. In the 
tight-binding model with spectrum fl}, it is given by 

G R (e, k N , k t , K) = J2 6 ^,K+bGo (e, k) . 



in Fig[2] also accounts for Umklapp processes, in which 
q z and q z i differ by b. 



where G R (s, k) = (is — + i/2r)~ , 6 = 2nn/d is the 
reciprocal lattice vector, and k = (k|| , fc z ) . The screened 
Coulomb potential is given by the RPA series (Fig. [2}. 
Because the in-plane motion is free, the in-plane momen- 
tum is conserved at the vertices. However, it is the out-of- 
plane quasi-momentum rather than the momentum that 
is conserved at the vertices. In addition to normal scat- 
tering, where the incoming and out-going bosonic mo- 
menta of the polarization bubble are the same, the series 



J 



FIG. 2: RPA series for electrons on a lattice. The incom- 
ing and outgoing bosonic momenta can differ by an arbitrary 
reciprocal lattice vector. 



Summing up the geometric series, we get for the 
Fourier transform of the dynamic screened potential 



D H (u,q\\,q z ,q z ) = D (q\\,q z )S(q z - q' z ) 



Do(<l\\,qz)'tt R (u,q L )J2 b Do(ci\\,qz +b)S(q' z - q z 



l + E b D (q lb q z +b)n^(to,q) 



(32) 



Here. 



4-7re 2 

A>(q||,&) = 3—2 

yii ' Hz 



(33) 



is the bare Coulomb potential and H R (ui,q) is the ana- 
lytic continuation of the Matsubara polarization bubble 



/d k 
- — 3G0 (ie + iuj, k + q) G (is, k) . 
(2tt) 

(34) 

I 



where Go(«e,k) = (ie — + isgne/2-r) -1 and 
q = (qib Qz) • The potential at the position of the 
most efficient resonant center, z; = d/2, is obtained from 
Eq.(|32| as 



D R (u,, qi] ,d/2,d/2) =J^J d ^e^-^D R {^,q z ,q z ,) = ^ + j 



dq z D Il R A 1 
2tt 1 + U R A 2 ' 



r 



where A\ and A2 are the lattice sums which can be per- formed explicitly 
I 



v- * ~ / ,n 2-Ke 2 d 

Ai = 2 - ,Do(q||,«* + 6) : 

b 

A 2 = Y, D o(^Qz + b)e- M / 2 ^ 



smh(qnd) 



qu cosh(q||d) — cos(g z c?) 

4:ne 2 d cos(q z d/2) sinh(gi| d/2) 



cosh(<7|iG?) — cos(q z d) 



(35) 
(36) 



In a strongly layered metal, the polarization operator (|34[) depends on q z only weakly, because the interlayer 
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hopping J c is small and the electron dispersion is almost 
two-dimensional. If this dependence can be neglected 
completely, the integration over q z can be readily per- 
formed with the result 



D R (w,q||,d/2,d/2) = ^p-tanh (jy) coth 

where (complex) momentum k is defined by the following 
equation: 



cosh(fcd) = cosh(<7|id) 



2ire 2 



-n(w,<j,|)sinh(g|jd) . (38) 



If the effective range of the screened Coulomb potential 
is much larger than the lattice spacing, Umklapp scat- 
tering is strongly suppressed. Consequently, in the limit 
of K^d -C 1, where n\ — Ane 2 ^ is the (square of) screen- 
ing wave vector, and v% is the three-dimensional density 
of states, Eq. (f3"Tj) reduces to the continuum limit result. 
For finite values of K^d, the role of Umklapp processes is 
quite important. 

Eq. (|37[) can be reproduced in an alternative way. In the 
limit of J c — > 0, the motion of electrons between different 
layers is forbidden. Therefore, the problem is equivalent 
to the one of screening of an external charge by paral- 
lel conducting layers^. Although both approaches are 
equivalent for J c = 0, the RPA method is more general 
since, unlike electrostatics, it allows one to consider the 
case of finite J c as well. 

The polarization bubble can be calculated explicitly in 
the diffusive (lot -C 1) 



Il R (uj,q) = V 3 



D\\q 2 + 8Tj 2 sm 2 {q z d/2) 
Dug 2 + St J 2 sm 2 {q z d/2) - i 



(39) 



and ballistic (u> 3> r 1 ) limits 



Tl R (ui ) q)=i/3 1H , for VpQn 3> min { J c dq z , oj} . 

VFq\\] 

(40) 

Here Dm = v 2 f t/2 is the in-plane diffusion coefficient. 
The additional assumption of large q\ | , employed in 
Eq.(|3Q"]). will be justified later. 

The local spectral weight <S^~ e (u>) for the electron- 
electron interaction has two distinct forms of the lu de- 
pendence. Assuming that the motion of electrons is bal- 
listic and strictly two-dimensional, we expand the imag- 
inary part of the potential in Eq. (|37p to first order in 
u/vFq\\ and perform the integration over gn. Typical 
value of qu turns out to be large: of order of the 2D 
screening wave vector K2 = 27re 2 ^2, where V2 is the 2D 
density of states. On the other hand, typical value of 
lu are determined by the temperature. Therefore, the 
assumption of lu /vpq\\ is satisfied for all reasonable tem- 
peratures: T <C vf&2- As the result comes from the first- 
order term in cj/vpqu, the resulting spectral density is 
linear in to 



S^iuj) = -nsu, (41) 
where t\b is the effective coupling for the electron-electron 
interaction 



>1B 



VF 



-9B (K2d) 



(42) 



and 



9B [X) 



1 coth 2 
dy 



coshy ( y + a; tanhf) 



1/2 



s coth |) 



3/2 



2.T 



for i<1 
for x ^s> 1. 



(43) 



r 



Large in-plane momentum transfers also help to jus- 
tify the assumption of two-dimensionality. Indeed, when 
calculating the polarization operator, one has to com- 
pare the difference of dispersions for the 2D and c-axis 
motions, i.e., 8e\\ = e k|| + q| | - e k| | ~ v F q\\ ~ v f k 2 
and 5e z = €k z +q z ~ £q z ~ Jdq z . Typical values of q z 
are of order minl^,^ 1 }- Hence, neglecting the z- 
component of the electron dispersion is justified as long 



as J <C vf^2/ 'dmin {k2, d 1 } , which is the case for any 
real layered metal. In the CL terminology 43 , Eq. flTj) 
corresponds to Ohmic regime. Substituting Eq. (|4ip into 
Eq. (j24|) and cutting off the ultraviolet logarithmic singu- 
larity in the lu integral at some high-energy scale Eq , we 
evaluate the integral for large E t as 



J 



'"° duj 



(1 — cos(wi)) coth — i sm(wt) 



= ]n(E Q t) - — sgn(t) + 7 + 



(44) 
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where 7 = 0.577... is the Euler's constant. Rescaling t 
by T -1 in the remaining integral of Eq. (|24]) . we find that 
the Ohmic regime of the spectral function corresponds to 
a power-law scaling of the resonant conductivity 

(t„(T)/<t A = C(t, b )(^\ \ (45) 

where the regularization dependent prefactor C (x) (for 
a hard cutoff regularization of the frequency integral) is 
given by 

C{x) = cos (^) e-T* . (46) 

As the distance between the planes increases, the tun- 
neling exponent ijb [cf. Eqs. (|42|43p ] decreases from 
its nominal value, given by the dimensionless electron- 
electron coupling constant in a bulk metal, to zero, which 
is to be expected. 

The expansion of the polarization bubble in w/vfQ\\i 
which led us to Eq. (|4"S"|) works as long as the subsequent 
integral over qu converges (and typical values of g|| are 
determined by the ultraviolet parameters of the theory) . 
This is the case for the leading term in the expansion; 
however, the next order term diverges logarithmically in 
the infrared. In the ballistic limit, the divergence results 
in a subleading cj 2 ln|w| correction to the spectral func- 
tion, which is not relevant at low temperatures. However, 
one of the factors of uu is replaced by r _1 in the diffusive 
limit (ui <C t -1 ), and the subleading term becomes larger 
than the leading, linear-in-w term. The resulting win |o;| 
behavior of the spectral function can be obtained accu- 
rately by starting with the diffusive rather than ballistic 
form of the polarization bubble [Eq. ([39|) ]. In the diffusive 
regime, one can still neglect the g z -dependent terms in 
11^ (w, q) , which amounts to an assumption of purely 2D 
diffusion, as long as lo 3> u>i, where 

^ = jTK3d {l, K2 d»l. ^ 

The resulting spectral weight is a sublinear function of 
the frequency 

S- c H^ DW logQ£), (48) 

where 

„ n - \ J ** d l 2 , n 2 d « 1 , . 

^ = {n 2 d D]] /P, K2 d»l, ( 5 °) 

and I = vpT is the mean free path. The corresponding 
conductivity 

a rcs (T) = a cl exp (tj 2D log (Tr) log (^f) 1 (51) 



increases with temperature faster than any power-law. 
Interestingly, the temperature dependence of the res- 
onant conductivity is similar to that of the zero-bias 
anomaly (ZBA) in the 2D diffusive case^ML. The dif- 
ference between the two results is in the dimensionless 
parameters: the dimensionless conductance, g — v 2 D\\ 1 
that controls ZBA, is replaced by 7720 for the resonant 
tunneling case. 

For frequencies smaller than u)\, cf. Eq.[57], typical 
gil's are large, which means the diffusion approximation 
breaks down again and the ballistic one should be used 
instead. We thus conclude that the resonant tunneling 
conductivity is given by the diffusive limit [Eq. (|51[) ] for 
temperatures in the interval uj± <C T <C 1/r and by the 
ballistic limit [Eq.([45p] for the rest of the temperatures. 

III. DISCUSSION OF THE RESULTS 

Several comments are in order. 

i) Formally speaking, the power-law scaling of the 
electron-assisted tunneling conductivity in the ballistic 
regime [Eg. (|4"5])] saturates for T > E . The value at sat- 
uration is the same as for the electron-phonon case: cr G i- 
However, Eq is of order of the plasma frequency, so that 
the electron-assisted mechanism of tunneling leads to a 
growth of conductivity for all reasonable temperatures. 
This feature may be used in experiment to separate the 
electron- and phonon-assisted mechanism: because the 
saturation temperature for the phonon mechanism may 
be not too high, the growth of the total conductivity up 
to the highest temperatures is indicative of the electron- 
assisted mechanism. 

ii) In contrast to the case of the interaction correc- 
tions to tunneling and transport conductivities in the 
ballistic regim o 50 i 51 , which are determined by the inter- 
action on a large spatial scale (of order of the ballistic 
Thoulcss length Vf/T), the electron-assisted tunneling 
conductivity is determined by the interaction at small 
distances (of order of the screening radius). At these dis- 
tances, the RPA works only for weak interactions, i.e., for 
e 2 /vF <§C 1, and Eq. (|4"B"|) is valid, strictly speaking, only 
in the perturbative regime. It is reasonable to assume, 
however, that if the interaction is not weak, the tunnel- 
ing exponent in Eq. (|45p is replaced by a non- universal 
quantity of order unity. 

iii) In a real system, both assisted mechanisms-the 
electron-phonon and electron-electron ones-operate si- 
multaneously. As we have explained before, the electron- 
phonon interaction is strong while the electron-electron 
one must be considered to be weak. To the extent that 
one can neglect the mutual influence of these two inter- 
actions, the total spectral function in Eq. ([M|) is a sum 
of two contributions S m (w) = S^p 11 + S°~ c (uj). Since 
Vb *C 1, one can expand the total conductivity (T*°| to 
first order in S!^ c (lu), which yields 

a^ = a ie Jl + T, B ln^-), (52) 
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where cr rcs is the electron-phonon contribution given by 
Eq. (f3"Tj) . Ignoring the logarithmic temperature depen- 
dence and using the low-temperature result for er ros [first 
line in Eq. (j3"Tj) ]. one sees that the electron-electron contri- 
bution dominates over the electron phonon-one for tem- 
peratures below a characteristic temperature 

T L = VW-Vd ■ (53) 

Using Eq. pTj) for A, estimating t]b as e 2 /vF and also 
the deformation-potential constant as A ~ Ms 2 ~ e 2 /ao, 
where M is the atomic mass and ao is the lattice con- 
stant, we obtain Tl ~ V 's/vf^d ~ 10 K. For tempera- 
tures above Tl but below the saturation temperature for 
the electron-phonon mechanism T s = \ujl>, the electron- 
phonon mechanism dominates. For T >T S the electron- 
electron interaction wins over again. Once the electron- 
electron mechanism becomes the dominant one, the log- 
arithms following the lowest order one in Ea.(|52p can be 
summed up into a power law. Therefore, a signature for 
the electron-electron mechanism is a power-law scaling of 
the conductivity (with an exponent of order unity) both 
at low and high temperatures, combined with an absence 
of saturation at high temperatures, as it was explained 
in item i) of this Section. The power-law increase of the 
conductivity with an exponent close to one has been ob- 
served in graphite both at low and high temperatures^. 

iv) Another mechanism, competing with those consid- 
ered in this paper, is the zero-bias anomaly (ZBA), i.e., 
the interaction corrections to the tunneling density of 
states. This mechanism leads to the temperature depen- 
dence of the resonant tunneling conductivity via the tem- 
perature dependence of the tunneling widths even in the 
absence of assisted processes. The interaction correction 
to the tunneling density of states for a 2D electron system 
can be divided into parts. The first contribution, 8v c i (e) , 
is determined by the interaction of electrons in the ab- 
sence of impuritie s 52 ' 53 . In a layered metal, 8v c \ (e) is 
energy-independent for e below the plasmon gap 53 , which 
can be safely assumed to be the case, so Sv c \ (e) does 
not contribute to the T-dependence of the tunneling con- 
ductivity. The second contribution, S1/4 (e) , comes from 
the interplay of electron-electron and electron impurity 
scatterings. In the perturbation theory, 5vd (e) is propor- 
tional to the inverse dimensionless conductance, (Eft)~ 
times the logarithmic function of g 33 i 50 i 53 both in the 
ballistic and diffusive limits. On the other hand, the 
electron-assisted tunneling conductivity in the ballistic 
regime [Eq. ([i5)) ] allows for an expansion in tjb^Eo/T. 
Therefore, the electron-assisted mechanism wins over the 
ZBA one if r/s is not too small: t\b ^> (Eft) 1 . In the 
diffusive regime, the electron-assisted and ZBA mecha- 
nisms are, generally speaking, of the same order. 

v) In Ref. [32l |. the insulating-like temperature depen- 
dence of p c (T) in high-T c cuprates was ascribed to reso- 
nant tunneling through a single resonant level, at energy 
Ae away the Fermi level. No assisted tunneling mech- 
anisms were invoked: the temperature dependence was 
coming entirely from the thermal distribution of conduc- 



tion electrons. In this model, the T dependence of the 
tunneling conductivity is insulating-like for T < Ae and 
metallic-like for T > Ae. However, averaging over the 
energy levels eliminates the T-dependence. We do not 
think that an assumption of a single energy level is a very 
realistic one and therefore involve boson-assisted mech- 
anisms, in which the T-dependence survives even after 
averaging over the energy levels (see, however, Sec. IIVI) . 



IV. COMPARISON WITH EXPERIMENT 

If the total conductivity of a layered metal is a sum of 
the band and resonant-tunneling contributions, as speci- 
fied in Eq.Q, then, depending on the parameters of the 
Boltzmann and tunneling parts of the conductivity, the 
total resistivity may exhibit a variety of T dependences: 
purely metallic (for weak tunneling), purely insulating 
(for strong tunneling), minimum at low T, maximum 
at high T, and both minimum and maximum. Figure 
[3] shows some of these behavior for typical band and 
resonant conductivities. As an example, we consider a 
model with a band resistivity pb{T) = po + aT 3 , with 
po = 100 and a = 10~ 3 in arbitrary units. The pa- 
rameters for the phonon-assisted resonant conductivity: 
A = 10, ojd = 2QQK. We plot a total resistivity for dif- 
ferent values of the elastic resonant conductivity, a e \. 

The appearance of the maximum in p c (T) has already 
been explained in Sec. H) one naturally has a mini- 
mum in the conductivity (and a maximum in resistivity) 
when adding up increasing and decreasing functions of 
the temperature. The minimum in p c (T) arises if a Tes is 
larger than itb at low temperatures. If quantum inter- 
ference effects can be ignored, the low-T dependence of 
us is due to electron-electron interactions. In a Fermi 
liquid, <jb = <7% (l — ci\jT 2 t/Ef) , where Oi is the resid- 
ual conductivity due to impurities and the dimension- 
less coefficient a\j parameterizes the strength of Umk- 
lapp scattering. If the resonant-tunneling conductivity 
decreases slower than T 2 , the net resistivity shows an 
insulating behavior. In the electron-assisted mechanism, 
the insulating upturn must occur at low enough tempera- 
tures, if exponent t\b of the electron-assisted mechanism 
[Ea.(|45p] is less than two. The electron-phonon inter- 
action is a marginal case because the low-temperature 
exponent equals precisely to two [see Eq. (TT5)) , first line]. 
In this case, whether the insulating upturn occurs or not 
depends on the magnitude of the electron-phonon and 
Umklapp interaction, and also on the amount of disorder 
in the sample. 

Experiment shows a variety of behaviors in p c (T). 
For example, p c i) is purely metallic in over doped 
cuprates; ii) has a minimum in under doped ones*; 
iii) is purely insulating in TaS^; iv) has both a min- 
imum and a maximum in graphite—; and v) has a 
maximum in Sr 2 RuOA K-(BEDT-TTF) 2 -Cu(SCN)22i, 
(Bio. 5 Pbo.5)2Ba 3 Co20 y and NaCo204 5 -. In our opinion, 
the most remarkable behavior is the one with a maximum 
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FIG. 3: resistivity in arbitrary units vs. tempera- 
ture for various values of resonant conductivity [a e \ = 
0, 0.00175, 0.003, 0.2). On the onset a non-monotonic low 
temperatures dependence for cr i = 0.2. 



in p c (T) . Whereas the insulating upturns can, in prin- 
ciple, be ascribed to phase transitions, which open gaps 
over the parts of the Fermi surface, the metallic behav- 
ior of p c (T) at low temperatures shows unambiguously 
that we are dealing with a well-defined metallic state. 
In all cases cited above, the insulating behavior at higher 
temperatures is not associated with, e.g., a ferromagnetic 
transition, as it is the case in manganites^. A mecha- 
nism explaining such a behavior without invoking metal- 
insulator transitions is suggested in this paper. In what 
follows, we focus on two of the materials with a max- 
imum in Pc - Sr 2 Ru0 4 and k-(BEDT-TTF) 2 Cu(SCN) 2 - 
and show that the data for these compounds can be fitted 
with our model 

As we have shown above, both the phonon-and 
electron-assisted mechanisms mechanism lead to an in- 
crease of the tunneling conductivity; the differences be- 
come important either at low enough or high enough 
temperatures. Given the uncertainty in other param- 
eters of the model, we perform the fit only for the 
electron-phonon mechanism. The low-temperature re- 
sistivity is dominated by ordinary Boltzmann transport 
and, in principle, may be calculated microscopically for 
a given Fermi surface. However, a large number of un- 
known quantities, such matrix elements for scattering 
and a complicated band structure, make such an ap- 
proach impractical. Instead, we extract the band part of 
the resistivity, pf , from the low-temperature part of the 
experimental data. We use the intervals of strictly metal- 
lic T-dependences between 10 and 50 K for Sr 2 RuC>4 and 
between 40 and 75 K for k-(BEDT-TTF) 2 Cu(SCN) 2 ), 
and then extrapolate the extracted dependence of p^ 3 to 
higher temperatures. The resonant part of the conduc- 
tivity is calculated numerically using Eg. (f!?5]) . The fit to 
the data for Sr 2 Ru0 4 and k-(BEDT-TTF) 2 Cu(SCN) 2 is 



40 




50 100 150 



T(K) 

FIG. 4: p c vs temperature. Solid: experimental data; dashed: 
fit into the phonon-assisted tunneling model. Top: Sr2Ru04. 
Fitting parameters: a c \ = 47 • 10 3 cm -1 , loo = 57 K and 
A = 18.5. Bottom: k-(BEDT-TTF) 2 -Cu(SCN) 2 . Fitting pa- 
rameters: Oc\ = 1.5 fi" 1 cm" 1 , ujd = 140 K, A = 18.9. 



shown in the top (bottom) panels of Fig. [4] The quality 
of the fit and reasonable values of the parameters suggest 
that the phonon-assisted model is a viable mechanism of 
the c-axis anomaly at least in these compounds. The data 
for Sr 2 RuC>4 shows a tendency to saturation for T > 400 
K, which is expected for the electron-phonon mechanism. 
As the electron-assisted mcchanism-not included in the 
fit-would lead to a further decrease in the resistivity, we 
can only speculate that effective coupling for this mecha- 
nism is very small in this material; this can be related to 
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a rather large distance between the planes. At tempera- 
tures above 700 K (not shown in Fig. rjj, the resistivity 
starts to rise slowly again. Although a re-entrant metal- 
lic behavior is not explained directly by our model, it 
can be understood if one recalls that for temperatures 
above the bandwidth of the resonant levels, the impu- 
rity band can be effectively replaced by a single level. 
As we have explained in Sec lIIIi in a single-level model 
the thermally activated resonant-tunneling conductivity 
is insulating-like for temperatures below the energy dif- 
ference between the resonant level and metallic-like for 
higher temperatures. 



V. SUMMARY 

In this work, we suggested an explanation of the 
non-metallic temperature dependence of resistivity, ob- 
served in various layered metals. This explanation is 
based on the interplay between two conduction chan- 
nels of transport: the band one with metallic-like tem- 
perature dependence and the resonant-tunneling one 
with the insulating-like temperature dependence. We 
developed a theory of electron-assisted tunneling that 
complements the previously known theory of phonon- 
assisted tunneling. According to our picture, the low- 
temperature part of the resistivity is determined by re- 
laxation of quasi-momentum due to the electron-electron 
and electron-phonon interaction. At higher tempera- 
tures, the electron-assisted and phonon-assisted mecha- 
nisms increase the probability of resonant transmission 
and lead to a decrease of resistivity with temperature. 

Our model relies heavily on the assumption that reso- 
nant sites are indeed present in real materials. Although 
there is an evidence that inter-plane disorder does en- 
hance the c-axis conductivity^!, a more direct verification 
of the resonant-level hypothesis is needed at the moment. 
A combination of experimental techniques, in which dis- 
order is introduced controllably and detected spectro- 
scopically, with first-principle computational techniques 
should help to resolve the issue. Such an approach which 
combines low-dosage intercalation of graphite with first- 
principle calculations of the energy levels of intercalated 
impurities is currently being pursue d 48 ; 55 . 
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APPENDIX A: DERIVATION OF EQJ3TI 



be extended ro infinity 



/(*) 



dw- 



o 

" D diO 
~ 2~ 



T 

D 



(1 - cos(wt)) ( 



coth 1 

2T 



ui 



w (l - e lut ) 



(Al) 



D 



For T < u D /V\, one can expand the exponential in 
Eq.© as follows 



e- x ' 2 / dt 



Using identities 

x . i 



1 + A 



duiuje 1 ' 



sinh (irx + iO) 
xe lpx dx 



-dx 



and 



7T 



sinh {ttx + iO) tt 



= *(coth(f 



1 sinh 



(A2) 

-2 (V 



we get 
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Performing integration over frequency, we derive the low- 
temperature asymptotics of a les 
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Next, we consider the case of ujd/^/X <c T -C ujd- Per- 
forming integration over frequency, we find an explicit 
formula 
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U! D t 



I ir 2 T 2 
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TT 2 T 2 



D uo' z D t' A LU 2 D smh z (irTt) 



(A3) 



The time integration in Eq. (|29p can be performed in the 
saddle point approximation, where the saddle point so- 
lution is to be found from 



Xf'(t) - 27rTcoth(7rTt) + - = . 



(A4) 



Recalling that T<wd, Eq.JA4j) can be simplified further 
2cot(w) uj 2 2 (ud\ 2 , s /. s 



where y = —iirTt. Solving Eq. (|A5j) to leading order in 
u>d/T, we find 

i 

t* = — . 
2T 

As a result, we find that for ujd/VX <C T <C ujd the 
resonant tunneling conductivity is given by 



At low temperatures T « wj), the upper limit of the 
integration over frequency in second term of Eq. (|28p can 
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For high temperatures (T ^> u>d), the function fit) For Xurj/T <C 1, the exponent in Eq. (|A7|) can be ex- 
can be approximated by panded leading to 

/(f)~^(**T-it). (A6) 

The resonant conductivity in this temperature range is 

givenby s= -/ a 

cr e i 7 sinh 2 (vri + 2 0) \ 3T v 7 = i _ ^£ ( A 9) 

(A7) 9T ' v ; 

For \ujd/T ^> 1, the integral in Eq. (|A7|) is evaluated by 
the saddle point approximation, yielding 



Crcs / 37tT ^ Awd^ 



4A^ eX H^J- (A8) 



1 S. L. Cooper and K. E. Gray, in Physical Properties of High 
Temperature Superconductors, edited by D. M. Ginsberg, 
(World Scientific, Singapore, 1994), p. 61. 

2 I. Terasaki, Y. Sasago, and K. Uchinokura, Phys. Rev. B 
56, R12685 (1997). 

3 S. M. Loureiro, D. P. Young, R. J. Cava, R. Jin, Y. Liu, P. 
Bordet, Y. Qin, H. Zandbergen, M. Godinho, M. Nuiiez- 
Regueiro, and B. Batlogg, Phys. Rev. B 63, 094109 (2001). 

4 I. Tsukada T. Yamamoto, M. Takagi, T. Tsubone, S. 
Konno, K. Uc hinokura J Phys Soc Jpn. 70, 834 (2001), 
|cond-matt/0012395| 

5 T. Valla, P. D. Johnson, Z. Yusof, B. Wells, Q. Li, S. M. 
Loureiro, R. J. Cava, M. Mikami, Y. Mori, M. Yoshimura, 
T. Sasaki, Nature 417, 627 (2002). 

6 Y. Maeno, H. Hashimoto, K. Yoshida, S. Nishizaki, T. Fu- 
jita, J. G. Bednorz, and F. Lichtenberg, Nature 372, 532 
(1994). 

7 N. E. Hussey et al., Phys. Rev. B 57, 5505 (1998). 

8 A. W. Tyler, A. P. Mackenzie, S. NishiZaki, and Y. Maeno, 
Phys. Rev. B 58, R10107 (1998). 

9 W. J. Wattamaniuk, J. P. Tidman, and R. F. Frindt, Phys. 
Rev. Lett. 35, 62 (1975). 

10 see N. B. Brandt, S. M. Chudinov, and Ya. G. Ponomarev, 
Semimetals: I. Graphite and its compounds, (North- 
Holland, Amsterdam, 1988) and references therein. 

11 J. Singleton and C. Mielke, Contemp. Phys. 43, 63 (2002). 

12 N. Kumar and A. M. Jayannavar, Phys. Rev. B 45, 5001 
(1992). 

13 P. W. Anderson, The theory of superconductivity in the 
high T c cuprates (Princetion University Press, 1997), p. 
50. 

14 N. Dupuis, Phys. Rev B 56, 9377 (1997). 

15 Abrikosov WL 

16 P. Wolfle and R. N. Bhatt, Electron localization in 
anisotropic systems, Phys. Rev. B 30, 3542 (1984). 

17 D. B. Gutman an d D. L. Maslov, Phys. Rev. Lett, (in 
press): larXiv:0704.0613l 

18 R. E. Prange and L. P. Kadanoff, Phys. Rev. 134, A566 
(1964). 

19 J. Singleton, P. A. Goddard, A. Ardavan, A. I. Coldea, S. 
J. Blundell, R. D. McDonald, S. Tozer, and J. A. Schlueter, 
Phys. Rev. Lett. 99, 027004 (2007). 

20 L.B. Ioffe, A.I. Larkin, A. A. Varlamov, and L. Yu, Phys. 



Rev. B 47, 8936 (1993). 

21 P. Moses and R. H. McKenzie, Phys. Rev. B 60, 7998 
(1999). 

22 A. Millis, Nature 417 599 (2002). 

23 Xu Du, S.-W. Tsai, D. L. Maslov, and A. F. Hebard, Phys. 
Rev. Lett. 94, 166601 (2005). 

24 Experiment 2 ^ shows that the temperature dependence of 
the transport scattering rate in graphite obeys l/r tr = 
0.065T for 25 K< T < 200 K. An excellent agreement of 
this result with the theory of the electron-phonon interac- 
tion indicates that the scattering mechanism is quasielastic 
scattering on phonons, when the quasiparticle lifetime, r 9 , 
and transport time are the same. At T = Tm = 40 K, 
we then have TuT q ~ 15. Notice also that since scattering 
is quasielastic, the Boltzmann equation would not break 
down for graphite even if the product Tmt, were less than 
one. 

25 U. Lundin and R. H. McKenzie,Phys. Rev. B 68, 
081 101 (R) (2003). 

26 A. F. Ho and A. J. Schofield, Phys. Rev. B 71, 045101 
(2005) 

27 J. Paglione, C. Lupien, W.A. MacFarlane, J.M. Perz, L. 
Taillefer, Z.Q. Mao, and Y. Maeno, Phys. Rev. B 65, 
220506(R). 

28 M. Turlakov and A. J. Leggett, Phys. Rev. B 63, 064518 
(2001). 

29 M. J. Graf, M. Palumbo, D. Rainer, and J. A. Sauls, Phys. 
Rev. B 52, 10588 (1995). 

30 A.G. Rojo and K. Levin, Phys. Rev. B 48, 16861 (1993). 

31 C. Uher, R. L. Hockey, and E. Ben- Jacob, Phys. Rev. B 
35, 4483 (1987). 

32 A. A. Abrikosov, Physica C 317-318, 154 (1999). 

33 B. L. Altshuler and A. G. Aronov, in Electron- Electron 
Interaction In Disordered Systems, ed. by A.L. Efros and 
M. Pollak (Elsevier, 1985), p.l. 

34 J. G. Analytis, A. Ardavan, S.J. Blundell, R. L. Owen, E. 
F. Garman, C. Jeynes and B. J. Powell, Phys. Rev. Lett. 
96, 177002 (2006). 

35 L.I. Glazman and R.I. Shekhter, Sov. Phys. JETP 61 163, 
(1988). 

36 N. S. Wingreen, K. W. Jacobsen, and J. W. Wilkins, Phys. 
Rev. Lett. 61, 1396 (1988); Phys. Rev. B 40, 11834 (1989). 

37 L.I. Glazman and M.E. Raikh, JETP Lett. 47, 452 (1988). 



14 



T. K. Ng and P.A. Lee, Phys. Rev. Lett. 61, 1768 (1988). 
A.I. Larkin and K. Matveeev, Sov. Phys. JETP 66, 580 

(1987) . 

A. K. Geim, P. C. Main, N. La Scala, L. Eaves, T. J. Foster, 
P. H. Beton, J. W. Sakai, F. W. Sheard, M. Henini, G. Hill, 
and M. A. Pate, Phys. Rev. Lett. 72 2061 (1994). 
L.I. Glazman and K.A. Matveev, Sov. Phys. JETP 67 1276 

(1988) . 

D. Ephron, M.R. Beasley, H. Bahlouli, K.A. Matveev 
Phy.Rev. B 49 2989 (1994); H. Bahlouli, K.A. Matveev, 
D. Ephron, M.R. Beasley Phy.Rev. B 49 14496 (1994). 
A.O. Caldeira, A. J. Leggett, Phys. Rev. Lett. 46, 211 
(1981). 

H. J. Lipkin, Quantum mechanics; new approaches to se- 
lected topics, (North-Holland, New York, 1973). 

I. G. Lang and Yu. A. Firsov, Sov. Phys. JETP 16, 1301 
(1961); Sov. Phys. Solid State 5, 2049 (1964). 

G. D. Mahan Many-particle physics (Plenum Press, New 



York, 1981). 

47 A. M. Finkel'stein, Sov. Phys. JETP 57, 97 (1983); ibid 
59, 212 (1984); Sov. Sci. Rev. A 14, 1 (1990). 

48 S. Tongay and A. F. Hebard, private communication. 

49 L. B. Ioffe and A. J. Millis, Phys. Rev. B 58, 11631 (1998); 
ibid. 61, 9077 (2000). 

50 A. M. Rudin, I. L. Aleiner, and L. I. Glazman, Phys. Rev. 
B 55, 9322 (1997). 

51 G. Zala, B. N. Narozhny, and I. L. Aleiner, Phys. Rev. B 
64, 214204 (2001) 

52 D. V. Khveshchenko and M. Reizer, Phys. Rev. B 57, 
R4245 (1998). 

53 E. G. Mishchenko, A. V. Andreev Phys. Rev. B 65, 235310 
(2002). 

54 M. B. Salamon and M. Jaime, Rev. Mod. Phys. 73, 583 
(2001). 

56 H.-P. Cheng, private communication. 



